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Abstract. The radial neoclassical fluxes of electrons in the 1/v— regime are 
calculated with relativistic effects taken into account and compared with those in 
the non-relativistic approach. The treatment is based on the relativistic drift-kinetic 
equation with the thermodynamic equilibrium given by the relativistic Maxwell- 
Juttncr distribution function. It is found that for the range of fusion temperatures, 
T e < 100 kcV, the relativistic effects produce a reduction of the radial fluxes which 
does not exceed 10%. This rather small effect is a consequence of the non-monotonic 
temperature dependence of the relativistic correction caused by two counteracting 
factors: a reduction of the contribution from the bulk and a significant broadening with 
the temperature growth of the energy range of electrons contributing to transport. 

The relativistic formulation for the radial fluxes given in this paper is expressed 
in terms of a set of relativistic thermodynamic forces which is not identical to the 
canonical set since it contains an additional relativistic correction term dependent 
on the temperature. At the same time, this formulation allows application of the 
non-relativistic solvers currently used for calculation of mono-energetic transport 
coefficients. 
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1. Introduction 

The role of relativistic effects in hot plasmas has been recognized as important not 
only in astrophysics [U|2] but also in fusion, in particular, for the population of highly 
energetic runaway electrons in tokamaks [3]. However, the relativistic effects do not 
necessarily require an extremely high temperatures since they can be non-negligible 
even if T e is only on the order of tens of keV, i.e. T e <C m e0 c 2 . These effects appear due 
to the macroscopic features of the relativistic thermodynamic equilibrium given by the 
Maxwell- Jiittner distribution function [HE] . An example of such effects provided by the 
Maxwell- Juttner distribution function is given in a recent paper [5] , where the stability 
criterion for collisional heat transfer from hot electrons to ions with respect to the 
Coulomb decoupling is studied and it is found that relativistic effects lead to qualitative 
changes in stability criteria. While in non-relativistic plasmas criterion is given by 
T e /Tj < 3, relativistic effects makes it temperature-dependent and for T e>i > 75 keV the 
collisional coupling between electrons and ions becomes absolutely stable. 

Relativistic effects in fusion are surely not important for the ions, but the transport 
physics for electrons needs to be examined carefully for fusion reactor projects such as 
ITER [6HE] and DEMO [9lfT0] . in which the expected electron temperature is sufficiently 
high, T e ~ 20 - 50 keV, and for future aneutronic fusion reactors with D- 3 He and may 
be p- u B reactions, which require temperatures of up to 70 - 100 keV pTTHT3] . However, 
all transport codes (see, for example [H]) developed to date and applied for simulations 
of reactor scenarios are based on the non-relativistic approach. Furthermore, there is no 
quantitative definition of an applicability range for the non-relativistic transport models 
so far. 

Relativistic kinetics and MHD in plasmas are usually treated in the covariant 
formulation [H[2]. For neoclassical transport, however, the covariant formulation is 
not necessary since Lorentz invariance is of minor importance with respect to the 
characteristic drift velocity, V^jc <C 1. For this purpose, one can directly apply the 
relativistic drift-kinetic equation [15] with the relativistic Coulomb operator [4]. 

In this paper, the relativistic effects in the radial fluxes in the l/V-regime, which 
might be the most dangerous regime for future burning plasmas in stellarators and 
where the radial electric field plays no significant role, are estimated. This case was 
chosen for investigation because the role of the highly energetic tail of the distribution 
function in transport processes in this regime is expected to be the largest in comparison 
with other regimes. Indeed, the diffusion coefficient in the 1 / u-regime scales roughly as 
Vfc/v e oc v 7 , while in the tokamak banana-regime it scales as (? ce v e oc v~ l (here, is 
the radial drift- velocity, p ce is the Larmor radius and v e is the collision frequency). 

In Sec. [21 the relativistic drift-kinetic equation (rDKE) in the mono-energetic 
approach with a set of thermodynamic forces which differs from the canonical one is 
formulated. Only radial gradients are taken into account while the parallel electric field 
is excluded from consideration. In Sec. [3], rigorous expressions for the radial electron 
fluxes and transport coefficients in the \jv— regime are derived. In particular, the 
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expression for the relativistic radial heat flux is obtained. As a guideline, the paper [Ttj] 
was used, where the same was calculated in the non-relativistic approach. In Sec. HI the 
numerical comparison of the relativistic and non-relativistic transport coefficients and 
radial fluxes is performed, and in Sec. |5]a brief discussion of the results is given. 



2. Mono-energetic drift-kinetic equation for relativistic electrons 

The electron radial fluxes in toroidal plasmas (except the Ware pinch) can be calculated 
from the relativistic drift-kinetic equation (rDKE) for the first-order distribution 
function f el in the mono-energetic approach [T7HT9] . Using on a magnetic surface, 
with flux-sufrace label p, the set of variables (s, u, A), where s is the coordinate along 
the field-line, u = p/m e0 = is the momentum per unit mass, 7 = ^/l + u 2 /c 2 is the 
Lorentz-factor, A = (1 — £ 2 )/fr is the normalized magnetic moment, where £ = u\\/u is 
the pitch and b = B/Bq is the normalized magnetic field with the reference field Bo, the 
mono-energetic rDKE can be written as 

V(U) - Mu)C(f el ) = - (v dr • v P ) (1) 

The first term in Eq. (CQ) is the mono-energetic Vlasov operator, V = (v\\h + Vd r ) • V s , 
where h = ~B/B and V s is the gradient within the magnetic surface (here, A = 0). 
The second term is the pitch-angle scattering operator with the deflection frequency 
z^d(iO = Vp(u) + Vd(u) (the complete expressions for relativistic v e £ and v e £ are given 
in 



Appendix A) and the Lorentz operator is 



C=U(*±). (2) 



b d\ V 9X 
The relativistic drift velocity can be written as 

^>' B - m C i '' B>?fi (3) 

with E = — V$ = — $'Vp and $' = d^/dp, where $ is the plasma potential (here 
and below, e = |e|). One can see that only the last term in Eq. ([3]) contributes to 
p = Vd r ■ Vp on the right-hand side (RHS) of Eq. ([T|). Since our treatment is limited 
to the l/i/-regime, only such values of E for which electrons with large E/vB make no 
significant contribution to transport are considered. In this case, the E x B drift term 
can be omitted in the Vlasov operator, i.e V — v»h • V s . (In the more general case, this 
term must be included to obtain the -y/zz-regime which is more complex for analytical 
treatment and is not considered here.) 

Thermodynamic equilibrium for relativistic electrons is given by the Jiittner 
distribution function [1] also known as the relativistic Maxwellian jl], which may be 
conveniently represented as 

f eM j(u,p) = - 1 ^ nf C M j(p r )e-^-V j (4) 

71 u te 
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where u te = pte/ m eo is the thermal momentum per unit mass with p te = y/2m^oT e and 
p r = m e0 c 2 /T e . The Maxwell- Jiittner distribution function is normalized by density, 
n e = J d 3 u feMJi and the normalization factor is 



C MJ (p r ) = .[^-^—- 1-^ + 0(1/ p 2 r ), /i r >l, (5) 
V 2p r A 2 (/i r ) 8p r 

where K n (x) is the modified Bessel function of n-th order. For convenience, the Maxwell- 
Jiittner distribution function is used in Eq. (JTJ with the Boltzmann-factor included: 



eMJ 



-e<S>/T e 



feMJ- (6) 



Since plasma parameters such as density and temperature only depend on the flux- 
surface label, p, the derivative in the right-hand-side of Eq. (CQ) can be expressed in 
terms of the thermodynamic forces, 

^l = [A 1 (p) + KA 2 (p)]F eMJ , (7) 

where k = p r (^ — 1) is the relativistic kinetic energy normalized by T e , and the 
thermodynamic forces A\ and A 2 are defined as 

MP ) = — - hj + ft Uf-— , (8a) 
^ 2 (P) = ^, (8b) 

J- e 

with ra' e = dn e /dp, T' e = dT e /dp, and the relativistic correction-term 

# 3 A 5 15 2 



^(/i r ) = p r [-±-lj--^ — + 0(1/ fl 2 r ), p r » 1. (9) 



Note that in contrast to the "canonical" set of the thermodynamic forces [TTHH], 
which depend only on the normalized gradients of density and temperature (n' e /n e and 
T e '/T e , respectively), and not on the absolute values of these plasma parameters, the 
first thermodynamic force A\(p) in the relativistic set Eq. (jSJ) contains an additional 
temperature-dependent term. 

Finally, the reduced mono-energetic rDKE can be represented as follows: 

(h ■ V s )/ el - ^EM. C(f el ) = -2- p [A^p) + kA 2 ( P )} F e Mj- (10) 

Note that similar to the non-relativistic formulation, the energy enters in Eq. (jTUl) 
only as a parameter in ^fV£>(u)/u and the solution of Eq. (ITOl) describes only the 
pitch- and spatial behavior of the distribution function f el , which is the same for both 
relativistic and non-relativistic approaches. With the proper choice of parameters and 
right-hand-side of Eq. ( fTUi) . the solution from such solvers as DKES [19] and NEO-2 [20] . 
which solve the non-relativistic DKE directly, can be interpreted as a solution of the 
mono-energetic relativistic DKE. 
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3. Relativistic radial fluxes 

In this chapter, the radial fluxes of particles and energy in the 1 / v— regime are calculated 
following Ref. [16] with the mono-energetic DKE treated in the relativistic approach. 

Equation ( flOl) can be solved by integration along the field-line. Here, only the 
trapped electrons, B /B m! ^ x < X < B /B m i n , are considered (B mSuX and B m[n are the 
absolute maximum and minimum of B on the given magnetic surface, respectively). 
Enumerating the local minima of B along the magnetic field-line by k and integrating 
Eq. ([TO]) over the bounce trajectory (assumed to be closed), one can obtain 

2 lMu) d ( Xl(k) df^\ = (k) dF eMJ 

u o\y ox J H dp K J 

with 

where 6pW is the radial displacement of an electron due to the magnetic drift after 
one bounce period. To solve Eq. (llljl . the following trick was used |16j . Applying the 
explicit expression for Vd r given by Eq. (EJ) to p — Vdr ■ Vp and using the fact that 
powers of £ = o~\/T — Xb with a = ±1 can be expressed as 

[m + 2)0 oX 

one can represent the integrand for 5p^ in Eq. ( fl2l) as follows: 

p_ \Vp\k G d f^ + e\ (14) 



£ 7 b 2 uj c0 d\ V 3 
where u c0 = eB /(m e0 c) is the cyclotron frequency, k G = n p • [h x (h • V)h] is the 
geodesic curvature of the magnetic field line and n p = Vp/|Vp| is the unit vector 
normal to the magnetic surface. Then 

= with fl*) = / ds ( (3 + a^2. (15) 

Using this relation and the fact that 1^ = = at the bottom of the magnetic 
wells (when A = B /B m i n ), the order of Eq. (fTTj) can be reduced, 

QfW H (k) u 2 



eMJ 



d\ 6A/( fc ) 7^d(w) dp 
The radial components of the particle and energy fluxes are given by 



(16) 



TP= (r e -Vp) = {^ j d 6 upf el ) , (17a) 

Q P e = (Qe • Vp) = I [ d 3 « e 0C 2 (7 - 1) P fel ) , (17b) 



where / e i is the solution of the relativistic drift-kinetic equation Eq. ( flQj) and (...) means 
averaging on the magnetic surface. 
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The conductive heat flux for relativistic electrons requires special attention. 
According to its physical definition [I31I8], the radial conductive heat flux can be found 
by extracting the advective and mechanical contributions from the radial component of 
the energy flux, 

<£ = <y e -T e rp-w e v, (18) 

where 

H-^Wt.-lJ/eM^g + K)^ (19) 

is the energy density related to the Maxwell-Jiittner distribution function, and V p = 
Y p e jn e is the flow velocity. Finally, the radial heat flux can be written as 



2 

Please note that this definition differs from the non-relativistic expression accepted in 
the neoclassical theory [T7t[T%] by the additional correction term 1Z. 

It is convenient to use common notations for both the particle and energy fluxes 
of the form Jj = (J hipf e id 3 u"), where J\ = T p and J 2 = Q p e /T e with h\ — 1 and 
/&2 = ac = p r {l — 1), respectively. Using in J d 3 u the variable k instead of u and 
performing the integration over A instead of pitch, J d£ = b/2J2 a f ^VICI) one can 
obtain from Eq. ( flTl) the following: 

r°° f _i_ i N 1 / 2 / r i/b (k ) ■ \ 

Then, substituting Eq. (JUj) into Eq. (1211) . an integration by parts over A can be 
performed. Finally, considering the averaging over the flux-surface as the limit of 
integration along the field-line and applying Eq. (JTBjl . the desired expression for the 
radial fluxes can be obtained, 

j rr r, e ~^ 5/2 / / 7+^ 5/ \ dhxF eMJ 

Ji = -n e G C M j dn — - — — — - I — - — hi , (22) 

Jo 1Vd{u) V 2 J dp 

with z>d(u) = V£,{u)/v e Q. One can check that the non-relativistic limit considered in 
Ref. [16] is recovered. 

The coefficient Go in Eq- ( |22l) . identical for both relativistic and non-relativistic 
formulations, contains all parameters for plasmas and magnetic configuration which are 
specific for the considered \ jv— regime, 

where R is the major radius and e e fr is the effective ripple amplitude (not shown here; 
for details see [H]). The expression for the radial fluxes in the \jv collisional regime 
calculated in the relativistic approach Eq. (1221) is, actually, the main result of this paper. 
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Figure 1. [Color online] a) Relativistic transport coefficients divided by 
corresponding non-relativistic values Lff are shown as the function of electron 
temperature. b) Integrand of L 2 2 (see Eq. (|25]l) is plotted for the different 
temperatures. The value T e = 1 eV is taken as the non-relativistic limit. 



Substituting in Eq. ( )22|) the derivative of the Maxwell- Juttner distribution function 
Eq. (CO), the radial fluxes can be expressed as 

Ji = -n e L H A i, (24) 

J=l,2 

where the thermodynamic forces A\ and A2 are defined in Eq. (jSJ) and the transport 
coefficients are equal to 

Uj = GqCmj / dK — — hihj, (25) 
Jo 1Vd{u) \ 2 J 

with i,j = 1,2. This definition satisfies Onsager symmetry. 

4. Comparison of relativistic and non-relativistic radial fluxes 

In this chapter, the role of relativistic effects in the radial neoclassical transport is 



examined. Using relativistic expression for ve>(u) (see Appendix A), direct numerical 
integration in Eq. (I25p can be done. For comparison, the expression for non-relativistic 
transport is appropriate. The latter can be obtained from Eq. ( 125]) by letting Cmj = 7 = 
1, k — v 2 /vf e , and with the non-relativistic expression for vu{v). Since the geometrical 
part of the transport coefficients is the same for both non-relativistic and relativistic 
approaches, the ratio of these quantities is a pure indicator of the relativistic effects. 

In Fig. Ufa), the ratios Lij/Lf* are shown as a functions of T e for i,j = 1,2 (here 
and below, the label "nr" indicates the non-relativistic quantities). One can see that 
the correction provided by the relativistic effects is not very strong (less than 7% for 
this range of temperature). However, a non-monotonic temperature dependence is not 
intuitively expected and requires an interpretation. 
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In Fig. QJb), the integrand for L 22 is plotted as a function of u/u te for different 
temperatures. One can see that a non-monotonic temperature dependence of L^jV^ 
can be explained by superposition of two counteracting relativistic effects. The first 
one appears due to a reduction of the contribution from the bulk of the distribution 
function and prevails in the low-temperature range, T e < 10 keV, leading to a decrease 
of transport coefficients (note that the slope in Fig. []7a) is almost the same for all 
transport coefficients in this temperature range). A decrease of the bulk contribution 
is caused by the specific feature of the Maxwell- Jiittner distribution function and can 
be estimated from Cmj Eq. ([5]). The second effect is caused by a broadening of the 
energy-range of contributing electrons and the shift of the maximum of the integrand 
into higher energies, and this leads to an increase of the transport coefficients with 
temperature. The latter effect appears to be important at higher temperatures. As one 
can see from Fig. H^b), in the non-relativistic limit (T e = 1 eV) the major contribution 
is coming from the electrons with u/u te ~ 1 — 4, while at higher temperatures this range 
becomes broader and for T e = 100 keV the corresponding range is u/u te ~ 1.5 — 6.5. 
This effect is weaker for Ln than for L 12 and L 22 due to the lower power of k in the 
integrand in Eq. ( 125]) . 

Unlike the non-relativistic case, the transport coefficients do not fully characterize 
the transport properties of a confined plasma (because of the relativistic factor 1Z in 
A\) and, consequently, the comparison of particle and energy fluxes is necessary as well. 
In order to make a comparison with the non-relativistic limit possible, let us consider 
two special cases: (a) n' e = $' = 0, and (b) T' e = 0, and the corresponding fluxes can be 
written, respectively, as 



J, 



(a) 



-Tl P 



J. 



(6) 



+ n Li i + l 



: 2 



t: 



n e L ul — - — 



(26a) 
(26b) 



In both cases, the ratio Ji/Jf T no longer contains the gradients and can be easily 
calculated. Note that in the case (6), Y p e /Y p e ,ra = Lu/L™, i.e. the relativistic correction 
for T p is identical to L n which is shown in Fig. Ufa). 

Following Eqs. (120]) and ( |26|) . the heat flux can also be represented in a similar 
manner: 



qP Xa)/ Te 



(4 + 211) L 12 + L 22 + 



q£>®/T e = -n e 



'12 



15 _ 
— + ATI + Tl 2 
4 



TL 



(27a) 
(27b) 



In Eq. ff27jl . the Onsager symmetry, L\ 2 = L 2 i, was used. 

In Fig. the ratios of r^/r^ nr , Q p JQ p e ' nv and q p /q p ' ni for both cases are shown. The 
same non-monotonic dependence as in the case of the transport coefficients is clearly 
indicated, and relativistic correction for the \jv radial fluxes is found to be less then 
10% for the temperature range checked. 




Figure 2. [Color online] The temperature dependence of the ratios Ji/J" r for two 
cases, (a) and (b), respectively, is shown. 



5. Summary 

In this paper, the neoclassical radial fluxes for hot electrons in the \ jv regime, which is 
specific to stellarators, has been calculated in the relativistic approach. The choice of 
ljv regime was motivated by a rather intuitive expectation that the role of relativistic 
effects in this regime should be most pronounced since the contribution to the radial 
transport from the tail of the distribution function is the largest (the integrand in the 
transport coefficients for the non-relativistic electrons in this regime scales in lowest 
order as oc v 7 ). 

For calculations, the reduced mono-energetic relativistic DKE was derived. Apart 
from the radial particle and energy fluxes, also the expression for the relativistic 
conductive heat flux was obtained. The definition for the radial fluxes in the relativistic 
approach has an important feature: the relativistic effects enter in the fluxes not 
only through the distribution function, but also through an additional temperature- 
dependent term in the first thermodynamic force. This relativistic term depends only 
on the temperature, in contrast to the canonical set of radial thermodynamic forces in 
which the logarithmic gradients of plasma parameters appear. Nevertheless, use of the 
proposed formulation has a big advantage: the transport coefficients with the relativistic 
effects taken into account can be calculated by the same numerical solvers which solve 
the non-relativistic DKE directly. 

Following Ref. [16], the radial fluxes were calculated from the relativistic mono- 
energetic DKE and the results obtained were compared with the corresponding non- 
relativistic quantities. It was found that the relativistic effects for hot electrons 
produce a modest, but systematic reduction of the radial transport (up to 10% within 
the temperature range relevant for fusion). However, a non-monotonic temperature 
dependence of the transport coefficients is somewhat surprising. This behavior is the 
result of two counteracting factors present for relativistic kinetics. The first factor is 
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related to a reduction in the relativistic Maxwellians of the weight of bulk electrons 
with an increase of the temperature. The second factor is caused by a broadening of 
the energy-range and a shift of the maximum contribution to higher energies. 

This initial investigation confirms the intuitive expectation of an absence of strong 
relativistic effects in the radial transport in stellarator fusion plasmas. At the same time, 
this conclusion is not general and a similar check must also be made for the banana- 
regime in tokamaks. Apart from this, maybe the most important task is the calculation 
of the parallel electron fluxes with the relativistic effects taken into account. Based 
on the results provided in this paper, one may expect that within the non-relativistic 
neoclassical treatment both the electron radial fluxes in the banana regime and the 
electron bootstrap current in hot plasmas are somewhat overestimated. 
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Appendix A. Relativistic expressions for deflection frequencies 



Expressing the deflection frequency for the test-particle a immersed in the background 
b through the diffusion coefficient for pitch-angle scattering, y^{u) = (2/u 2 )Z}^(m), and 
taking the general relativistic definition for Dgg(u) from Ref. [3], the general expression 
for Vd(u) with the relativistic Maxwellian can be written as follows: 

4 

Z/ff(tt) = V e0 C M j(Hr)—^X 
\/7T 



7 



it;' 



7 



7 



r /2 

7 



ii- 



r 



J0[2]02 



7 2 M 2^0[3]022 



f2 



II 



V 



(A.1) 



1 



c" u 

+ 



JO [2] 02 + 



8 c 2 . 

^T^j0[3]022 

7 Z u 



7' 



The specific functions ji[k]*(z) [3] are given by: 



(A.2) 



h[3m(z) = [-3*7 + (3 + 2z 2 )cr}/32z, 

where a(z) = ln(z + 7) with 7 = a/1 + z 2 and z = u/c. Since the leading order for 
z 1 is jo[2]o2 — z 2 /Q and jo[3]022 — z 4 /120, the non-relativistic limit, c — > 00, can be 
easily obtained [T71IT8] . 



1 

2^2 



erf(x) + 



erf'(x) 
2x 



(A.3) 



where x = v/v te 
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For i/p, the ion background can be taken as a non-relativistic Maxwellian. For 
calculations, it is sufficient to apply the high-speed-limit: 



One can see that the non-relativistic limit for this expression also exactly coincides with 
the classical approach. 
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